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in . Abstract 
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The phenomenon of wavepacket diffraction in space and time is investigated nu- 
Cd ' merically and analyticaUy, for a one-dimensional array of equally spaced finite-depth 

a; 

— ^ . wells. Theoretical predictions for the lattice at long times and at low scattering en- 

a: 

O ■ ergies, coincide exactly with the results for a single well. At intermediate and short 

. ^_, times compared to the classical passage time, the pattern shows both a broad diffrac- 

j^ ■ tive pattern and an interference pattern inside each diffractive peak. The diffractive 

structure persists for this case to infinite time. 
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1 Introduction 

In previous publications |jl], |^, ^ the phenomenon of wavepacket diffraction in space and time 
was described, both numerically in one and two dimensions [|I], 0|, as well as analytically 
for one and three dimensions at long times. ||^] The phenomenon was found in wavepacket 
potential scattering for the nonrelativistic, Schrodinger equation and for the relativistic Dirac 
equation. 

The essential characteristic of the phenomenon consists in the production of diffractive 
peak structure that travels in space, due to the interference between the incoming and spread- 
ing wavepacket and the scattered wave. The patterns are produced by a time independent 
potential. 

The peak structure exists for all packets, but, it survives only for packets that are initially 
narrower than a value related to the potential width (see below). For wider packets, the peak 
structure merges into a single peak. 

The effect is named wavepacket diffraction in space and time. 

For plane wave monochromatic waves, quantum scattering displays diffraction phenomena 
in time [Q induced by the sudden opening of a slit, or in space by fixed slits or gratings. 
The combined effect of time dependent opening of slits for plane monochromatic waves 
produces diffraction patterns in space and time.@ These patterns fade out asymptotically 
in time. Recent measurements of atomic wave diffraction0, have indeed demonstrated that 
the diffraction in time, process is supported by experiments. 

It was found in |l|, ^ |^ that there is no need to have a sudden (time-dependent) opening 



of a slit to generate a time-dependent diffractive train. The multiple peak traveling structure 
is generated by a time independent potential provided the packet width is narrower than the 
extent of the potential, either well or barrier. The time dependent behavior stems from the 
time-dependent spreading of the initial and scattered packets. 

The scattering events investigated so far considered only a single well. However, any 
target, be it a foil that is irradiated by some beam, or a quantum well or heterostructure 
scattering a current, involves a large array of wells. The purpose of the present work is to 
address the problem of wavepacket diffraction in the context of a one-dimensional array of 
equally spaced finite size and fine depth wells, as the time-honored Kronig- Penney model0. 
Our interest will be in dealing with a finite size array and not an infinite one. Therefore 
Bloch's theorem will not be of a help in treating the scattering event theoretically. The 
phenomenon we treat depends on the spatial extent of the wells. We will therefore avoid 
using 6 function type of potentials. It will be seen that, despite these complications, the 
outcome of the scattering is quite independent of the number of wells, their, depth, etc. 
The only aspect that will matter crucially is the width of a single well in the array. The 
results open up the possibility of using the phenomenon in solid state devices. In the next 
section we summarize some results of ref.j^]. Section 3 deals with the Kronig- Penney model 
analytically, and section 4 numerically. Section 5 brings a brief summary. 



2 Wavepacket diffraction in space and time 

1110 we found that, for t — i> oo it is possible to write analytical formulae for the diffraction 
pattern produced by a scattering event of a wave packet at a well or barrier. The formulae 
were derived for the simplest possible case that can be dealt with almost completely ana- 
lytically. This is the example of a wave packet scattering off a square well in one and three 
dimensions. The results, however, were shown to be independent of the shape of the well 
and the packet. 

As it is our intention to treat the one dimensional array of wells, we recapitulate here the 
one dimensional formula as a bridge to the lattice case. 

Consider a gaussian wave packet impinging from the left, on a well located around the 
origin. 



V{x) = -VoQ{w/2-\x\) (1) 

where, w is the width, is the Heaviside function, and Vq, the depth. Using the results 
of 0, the reflected and transmitted packets become 

^{x,t) = / (l){k,x,t) a{k,qo) e-' ^^^dk (2) 

where 0(a;, t, go) is the stationary solution to the square well scattering problem for each k 
and a{k, qo) is the Fourier transform amplitude for the initial wave function with average 
momentum Qq. 



As the effect is due to the interference between the incoming spread packet and the 
reffected wave we focus on the the backward direction for x < —w/2 



(p{k, X, t) = D{k, k') e' ^ (^-^") + F{k, k') e'' ^ (^+^'«+2 ^) 
D{k,k') = 1 

F(kk') - ^iM) 
^^^'^^ - Aik,k') 

E{k, k') = -2 i {k^ - k'^) sm{2 k'w) 

A{k, k') = {k + k'f e-2 '^' ^ -[k- k'f e^'^' "" (3) 

where 

aik, go) = e" ^ (^-^o)-<^' ('^-''o)' (4) 

with k! = Jk^ + 2 m |Vo|, and a, the width parameter of the packet. 

At fc ~ we have F k, —\. For long times and distances x » Xq » w we find using 
the properties of Gaussian integrals 



2 m TT __, , /m X , 9 x\ , 

2 W — - — e \sin{—-{xQ + 2i a' go) )| 



/2 m TT _^ / , m X xq . 2 a^ go m x, 

2 f^^ l-^ "I" -^oJ I 2 



^ = ^^ [ ^;, "^ + go^j (5) 

This expression represents a diffraction pattern that travels in time and persists. In H, 



we compared the expression above to the numerical solution. Excellent agrrement was found 
without resorting to any scale factor adjustment. (See figure 1 in 0). 
We further derived the condition for the pattern to persist to be 



a«J— (6) 

V % 

A relation between the initial width of the packet, the width of the well or barrier and 
the initial average momentum of the packet. Narrow packets only display the effect. 

We connected the results to the process of diffraction in tinieQ, and noted that the 
phenomenon of diffraction in space and time with wave packets did not need a time dependent 
opening of a slit, a time independent potential sufficed. It was found that the pattern did 
not fade out at long times as in the case of the diffraction in time|^, and it applied to any 
potential, be it repulsive or attractive. 

The key element for the obtention of the formula of eq.(^, was the replacement of F in 
eq.(|^), by the value at threshold, namely F = —1. This replacement is valid for long times 
at which the wild oscillations in the exponentials demand a very low momentum. Due to the 
fact that this value is independent of the type of well or barrier, the result holds in general 
for any kind of well or barrier. In all cases where the reflection coefficient F can be replaced 
by the value at threshold we obtain a receding multiple peak coherent wave train or, a single 
hump, depending on the initial width of the packet. This is the essence of the phenomenon 
of wavepacket diffraction in space and time. 



3 Kronig-Penney scattering 

Consider now, a one dimensional array of N square wells of width a to the right of the origin, 
a distance 5 + a apart from each other. We shift the location of the wells as compared to 
eq.(|^) for convenience. 

^-^ a a 

V{x) = -VoJ: e(--|x-/„--|) (7) 

where the left edge of the well is /„ = n (5 + a), while the right edge is located at 

' n t^ ~r tZ. 

In order to show that, this one dimensional Kronig-Penney lattice of finite size and depth 
wells yields essentially the same result as a single well, we will expand the reflection coefficient 
for the lattice to third order in k. 

In order to find the refiection coefficient, equivalent to F of eq.(^ we have to solve for 
the wave function as a function of k by means of plane waves, without resorting to Bloch's 
theorem that applies only for the infinite array. 

For each well we write 



^ = An e' '' "" + Bn e-' '' ^ r„_i < x < L 

= Cne"^^+ D^e-'^^ ln<x<rn 

= An+i e^ ^ ^ + S„+i e"^ ^ ^ r„ < X < Z„+i (8) 



with K = a/A;^ + 2 m Vq, m being the mass of the impinging particle. For the last well 
we demand, an outgoing wave, hence Bj^ = 0. 

Continuity of the wave function and its derivarive at the right and left edges of each well 
gives 
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where z\ = cos{k a) — i sin{K a) (1 + m^)/(2 u)^ Z2 = i sin{K a) (1 — 'U^)/(2 -u), u = K 



Due to the strong oscillations of e 



in eq.(0), the long time behavior of the wave 



function, is determined by the very low u = - regime. We therefore expand the matrices in 
powers of u. The reflextion coefficient F, will be then obtained from F = ^ as a polynomial 
in u. 

For a plane wave stationary solution with no imcoming waves from the rightmost edge 
of the wells and the number of wells A^ >> 2, we find after some lengthy algebra 
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where cq = cos{k,o a), X 



-(1 + 2 coA + (iV - 2) coA(l - e^ ' '') + 0{e') 



(10) 



-i sin{Ko a) (1 



o2 i k 5^ 



, e = — , and kq = a/2 m Vq is real 
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for wells, and imaginary for barriers. From the expressions above it is seen that to this order. 



the correction to F ^ — 1 originates from the term multiphed by (N-2). However this term 
is identical for both the expressions above. Therefore, to this order, Aq and —Bq are almost 
identical. 

We have carried out an expansion to O(e^) also. The expression becomes quite involved 
having some 20 terms each for Aq and Bq. Even to this order, the leading terms still give 
F^ -1. 

The Kronig-Penney lattice displays the same low k behavior of the reflection coefficient 
as a single well. 

In light of the results of the previous section, it follows that the Kronig-Penney lattice will 
generate scattered wave packets with the same characteristics as a single well. In particular, 
the phenomenon of diffraction of wavepackets in space and time will apply to the lattice 
also. In the next section we will show numerical calculations that support this claim, and 
also depict scattering events for different number of wells at finite times in contrast to the 
results found here for long times. 

4 Numerical results for the Kronig-Penney lattice 

The following series of pictures depict the amplitude of the wave function obtained with the 
numerical algorithm described in [jl|. We use Gaussian wells defined by 



Vix) = ~Voj:expU-^-^) (11) 

n=0 V / 



Where the equivalence to the square wells array of eq.(0) is given by 6 ~ ^. The 
impinging packet is a gaussian wave packet traveling from the left with an average speed f , 
initial location xq, mass m, wave number q = m v and width a, 

i! = C expU q {x - Xq) — - — j (12) 

We choose the parameters to be Vq = 1 a = 1, x„ = 8, m = 20, g = 1, a = 0.5, and 
xo = —20, that implement the constraint of eq.(^. In the case of a wider packet we find a 
smooth hump proceeding backwards and forwards P|. 

Figures 1,2,3 show, the backwards scattered, the forward scattered and the well region 
wave amplitudes at t=15000 for various numbers of wells. 

The figures show clearly that the diffraction pattern that dominates the wave profile 
in the backwards direction is the one of a single well, as expected from the results of the 
previous section. 

The amplitude of the wave is modulated by the diffraction pattern, with the interference 
pattern appearing as oscillations inside each broad peak. The number of oscillations is iden- 
tical to the number of wells, higher order interferences are blurred out into the background. 

In figure 4 we depict a detail of the wave inside the 12*/i well for the case of 49 wells as 
compared to the wave for the case of a single scatterer. The behavior is almost identical. 
The pattern repeats itself inside each well of the lattice, but the amplitude varies in space 
(well number), and time. 

In the previous section, we found that the long time behavior should be dominated by 
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Figure 1: Wave amplitude in the backward region for different numbers of wells. 
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Figure 2: Wave amplitude in the forward region for different numbers of wells. 
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Figure 3: Wave amplitude in the wells region for different numbers of wells. 
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Figure 4: Detail of the wave amplitude inside a certain well. 
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the diffraction pattern of a single well. Therefore, the interference ripples should diminish 
as time elapses. 

Figure 5 shows the results for a longer time and 49 wells. This calculation demanded al- 
most a day of computer work. Further increase in the time makes the calculation prohibitive. 
However, this increase should be enough in order to see if there is a trend of dissipation of 
the interference ripples. 

Indeed, the dissapearance of the broken peaks is most evident at long distance, at inter- 
mediate distances it is hard to assess if this is the case. We can expect that at even longer 
times the broken peaks will converge into the diffraction pattern humps broadening them to 
the size obtained with a single well. 

5 Brief summary 

We have found that the phenomenon of wavepacket diffraction in space and time remains 
intact for a one dimensional array of wells. On the experimental side, the present results 
suggest that besides the experiment simulated in [Q with Helium, it could be possible to use 
solid state devices, such as an array of quantum wells to implement the effect. Future work 
will focus on the more realistic situation of a higher dimensional lattice of scatterers. 
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Figure 5: Wave amplitude in the backward region for 49 wells at different times. 
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